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ABSTRACT 

We measure the clustering of Extremely Red Objects (EROs) in ps 8 deg 2 of the NOAO Deep Wide Field 
Survey Bootes field in order to establish robust links between ERO (z ~ 1.2) and local galaxy (z < 0.1) pop- 
ulations. Three different color selection criteria from the literature are analyzed to assess the consequences 
of using different criteria for selecting EROs. Specifically, our samples are (R-K s ) > 5.0 (28,724 galaxies), 
(I-K s ) > 4.0 (22,451 galaxies) and (/- [3.6]) > 5.0 (64,370 galaxies). Magnitude-limited samples show the 
correlation length (ro) to increase for more luminous EROs, implying a correlation with stellar mass. We can 
separate star-forming and passive ERO populations using the (K s - [24]) and ([3.6]- [24]) colors to K s = 18.4 
and [3.6] = 17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clus- 
tering properties and host dark halo masses, and cannot be simply understood as a single population. Based on 
the clustering, we find that bright passive EROs are the likely progenitors of > AL* elliptical galaxies. Bright 
EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the 
local Universe, making these the likely progenitors of > L* local ellipticals. This suggests that the progenitors 
of massive > AL* local ellipticals had stopped forming stars by z > 1 .2, but that the progenitors of less massive 
ellipticals (down to L*) can still show significant star formation at this epoch. 
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1. INTRODUCTION 

The pioneering / T -band near-infrared survey of 
lElston. Rieke. & Riekel (119881) was the first to identify 
galaxies with red colors of (R-K) > 5, and hypothesised 
these to be ei ther a new class of gala xy at z > 6, or ellipticals 
at z > 1- Graha m~& Devi (1 1996b named these galaxies 
Extremely Red Objects (EROs) for their extreme optical to 
infrared (IR) colors (e.g., (R-K) Vega > 5.0, (I-K) Vega > 4.0, 
(R - [3 .6])vega > 6.5) and recorded the first redshift for an 
ERO of z = 1 .44. We now know that the class of EROs 
consists of two z > 0.8 su b-populations: gala xies with old red 
stellar populations (e.g., Spin rad et al.l lT997) or star-forming 
galaxies with strong dust-obscuration (e.g., from interactio n 
and/or orientation; iDev et al.lll999HYan & Thompsonl2003T) . 

Since the highest mass galaxies at this redshift are of- 
ten found to be EROs (> 1 0" 5 M Q ; e.g.. iMoustakas et al.1 
I2004t IConselice et aLll2008l) . and bee ause EROs also exhibit 
very strong clustering, comparabl e to the local clustering 
of massive elliptical galaxies (e.g., iDaddi etal.ll2000l [2002; 
iBrown et al.l 120051) . EROs have been posited as the progen- 
itors of massive local ellipticals. This is important because 
the formation history of massive galaxies provides a crucial 
test of galaxy formation models. A key part of any such 
test is to link high and low redshift galaxy populations, and 
with the advent of large-area IR surveys, EROs provide an 
ideal and easily accessible testing ground for this purpose 
dGonzalez-Perez et al.ll2009h . 

One powerful means of Unking galaxy populations across 



redshift is through clu stering and the (dark matter) ha lo oc- 
cupation model (e.g., IMoustakas & Somervilld 12002b . In 
the case of massive dark matter halos, simulations of dark 
matter structure evo lution such as the Millennium Run 
(Sprin gel et al.l l2005) demonstrate that massive halos evolve 
rapidly in mass but little in (comoving) space over cosmic 
time. Since the most massiv e galaxies occupy the most mas- 
sive halos (e.g., Seliak 2000), they too grow in mass but their 
(comoving) spatial distribution evolves slowly. This means 
that accurate measurements of ERO spatial clustering and 
host halo mass can be used to link these massive z > 0.8 galax- 
ies with local galaxy populations. 
Previous studies have all found that EROs are very 



strongly clu s tered ( e.g., LD addi et al. 2000; Firth et al. 2002 
Roche et all l2002t iDaddiet alJ l2002t IBrown et aTT l2005t; 
Georgakak is et al.l 120051: iKong et all 120061 120091: iKim et al.1 
201 ll) . However, the existing estimates of the clustering 
strengths, as measured by the scale ro where the two-point 
correlation strength is unity, show a broad range (ro = 5.5- 
17/r'Mpc) that makes it difficult to link the higher redshift 
EROs with local galaxy populations. Many of these discrep- 
ancies are likely due to differing selection criteria and cosmic 
variance due to small survey areas and sample sizes. Differ- 
ences in methodology between studies, such as redshift dis- 
tribution models and uncertainty estimates, are also possible 
contributors. We discuss the consequences of these factors 
below. 

The diversity in ERO selection criteria are largely due to 
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pragmatic choices for selecting high-redshift galaxy popula- 
tions in different data sets. The present range of ERO selec- 
tion c riteria include (but are not limited to): (R- K)^^ > 5.0 



(e.g IBrown et alJ 120051) . {I-K) y e ^ > 4.0 (e.g.. IKong et all 
[2009b, (R ~ #)ve g a > 4.0 (e.g ., Firth etd] 12002b and (R^ 
[3.6])ve ga > 6.5 (e.g. lWilson et al.ll2004h . The consequences 
of using different ERO selection criteria have been considered 
in some studies (e.g.. lYan & T hompson 20 031 IWilson et all 
l2004t IConselice et al.l2008tlKim et aU201 1[) . IConselice et all 



(2008) found that ss 60% of EROs were selected by both 
(I—K) and (R-K) criteria, with the main differen ce bein g a 
shift in the mean redshift (f). Similarly, Ki m et all (1201 ll) at- 
tribute the difference in clustering between (r-K) and (i—K) 
selected EROs to b e a consequence of their difference in z. 
IWilson et all ( 120041) replaced the ZT-band with IRAC 3.6/xm, 
and suggested that ERO samples selected by this color will 
include more dusty EROs with star form ation than A'-band 
selected EROs (as we will discuss in $3.21 >. 

Morphological and spectroscopic studies of EROs show 
them to be a mixture of galaxies with and without star 
formation (e.g., ISmail et al.1 l2002t iMoustakas et al.1 l2004t 
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2006, IWilson et al.1 120071: IConselice et al.1 120081: 
2009). The relative fractions of dusty galaxies 



with active star formation (hereafter, star-forming), and galax- 
ies dominated by old stellar populations and little ongoing 
star formation (hereafter, passive) depends on the ER O selec- 
tion criteria and the limitin g magnitude (e.g., Moustak as et al.l 
l2004tlConselice et al.l2008|) . The sub-populations can be very 
broad, as star-forming EROs can have star formation rate s 
(SFRs) of between 10-200 M Q /yr (e.g. JMessias et al.l2010l) . 
but distinguishing further between sub-populations of star- 
forming EROs is beyond the scope of this study. Separat- 
ing passive from star-forming EROs is important because we 
should not expect the two ERO types to possess the same stel- 
lar mass, occupy the same mass dark matter halos or be the 
progenitors of the same local galaxy populations. 

Different fractions of these two ERO sub-populations se- 
lected by different ERO criteria may well explain the large 
differences seen in ERO clustering studies. Unfortunately, 
accurate morphological and spectroscopic classifications are 
often unavailable for the sources used in ERO clustering 
studies. In most cases, colors are used to distinguish be- 
tween s tar-forming and passive EROs. NIR colo rs (e.g., 
(J-K); Pozzetti & Mannucci 2000; Kim etlI] |20Tlh cannot 
adequately distinguish between the ERO types if they do not 
bracket or straddle the 4000A break, which does not work for 
all EROs (only z > 1 .4), a nd requires high precision photom- 
etry (e.g. JStern et al.ll2006l). T he mid-IR color te chniques in- 
troduced bv lSternet all (120061) ([3.6] - [24]) and IWilson et all 
(2007)([3.6]-[8.0]) are more promising because they use the 
emission from hot dust in star-forming galaxies to obtain a 
better separation of the two ERO types. However, these mid- 
IR color techniques require space observatory imaging, which 
is not always available. 

Accurate spatial clustering measurements require reliable 
redshift distributions, but th ere are relatively few EROs with 
spectroscopic redshifts (e .g., Grah am & Devil 1 996b iDey et al.1 
U999tlCimatti et al.ll2002T: IConselice et al.ll2008l) . and spectro- 
scopic samples at high-redshift are vastly incomplete (e.g., 
Coop eFet al.ll2007l) . Typically, redshift distributions are as- 
sumed to have a particular shape or are derived from photo- 
metric redshifts or both. This can be a problem because spa- 
tial clustering analyses are strongly dependent on the shape 



of the redshift distribution, contributing to the spread in pre- 
vious clustering estimates. Unfortunately, using photometric 
redshifts and models of the redshift distribution remains the 
best methodology available at present. 

An additional problem is that the angular correlation func- 
tion (w(9)) is usually fit as a fixed power-law, which may not 
correctly represent the da ta even if it is th e only option given 
the s ample sizes (e.g., iFirth et aT] 120021: iGeorgakak is et al.l 
120051) . Furthermore, the uncertainties in the correlation func- 
tion are frequently assumed to be Poisson with no correla- 
tions between data points (e.g.. IKong et al.ll2009t iKim et al.1 
2011). An incorrect power law, possibly combined with a 
poor uncertainy model can easily result in incorrect estimates 
of the spatial correlation length. A likely consequence of 
these problems is the broad range of ERO clustering estimates 
that are consistent with their having been the progenitors of 
anything from local L* ellipticals to A bell type clusters (see, 
iNorberg et al.l l2002: Aba di et al.ll 19981 for corresponding lo- 
cal galaxy population j*o values, respectively). 

In this paper we present a comprehensive measurement of 
ERO clustering. The sample sizes and survey area are large 
enough to address all these problems. We consider the largest 
ERO samples to date, selected in either the K s or IRAC [3.6] 
bands over 7.80 deg 2 or 8.53 deg 2 of the NOAO Deep Wid e 
Field Survey (NDWFS) Bootes field dJannuzi & Devil 19991) . 
respectively. The NDWFS Bootes field provides a factor of 
three larger survey area than any previous ERO clustering 
study. We study three ERO samples using the selection cri- 
teria, (R-K 5 ) > 5.0, (I-K s ) > 4.0 and (/- [3.6]) > 5.0, in an 
attempt to establish the differences and similarities between 
these samples and to better understand previous results. Our 
ERO samples contain 28,724 (R-K s ) > 5.0, 22,451 (I-K s ) > 
4.0 and 64, 370 (/- [3.6]) > 5.0 selected EROs, down to mag- 
nitude limits ofK s < 19.4, K s < 19.4 and [3.6] < 18.5, respec- 
tively. The sample sizes are 5 times larger than any previous 
study, a nd 10 times larger than our previous study of EROs in 
Bootes (IBrown et al.ll2005l) . 

We use optical NDWFS B^RI photometry dJannuzi & Devi 
1999), spectra from the Active Galactic Nuclei (AGN) and 
Galaxy Evolution Survey (AGES; Kochane k et al.ll2012l) and 
various small projects with instruments including Keck, Gem- 
ini and the Kitt Peak National Observatory Mayall 4m, 
in conjunction with IRAC photo metry from the Spitzer 
Deep Wide-Field Survey (SDWFS: lAshbv et al.l2009h . NEW- 
FIRM /fj-band photometry (Gonzalez et al. (in prep.)), 
and MIPS 24 fj,m photometry fro m the MIPS AGN and 
Galaxy Evolution Survey ( MAGES : IIannuzi et al.ll2010h . We 
test mid-infrared ERO separation techni ques utilizing equiv- 
alent B W RIK S a n d IRAC photometry ilBundy et all 120061: 
iDev et al.1 120071: iBarmbv et al.1 120081 MIPS photometry 
from the Far-Infrared Deep E xtragalactic survey (FIDEL; 
iDickinson & FIDEL team 20 07b and DEEP 2 spectroscopic 
indicators and ID spectra dYan et al.l 12009: Ne wman et al.l 
120 121 and references therein) in the Extended Groth Strip 
(EGS). 

This paper is structured as follows. We describe the data 
sets, photometry and photometric redshift estimation in $2] 
Our ERO samples are defined and our ERO separation tech- 
nique is described in ^3] In ^4] we present our clustering 
analysis and results, followed by a discussion and summary 
in ^5] Throughout this paper, unless otherwise specified, 
all magnitudes are in the Vega magnitude system. Spitzer 
IRAC and MIPS magnitudes are quoted using square brackets 
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(e.g., 3.6 fim is [3.6]). We introduce the following nomen- 
clature for our ERO samples: (R-K s ) > 5.0, (I-K s ) > 4.0 
and (7- [3.6]) > 5.0 selected ERO samples are respectively 
referred to as the EROra-, ERO/a- and ERO/3.6 samples. We 
assume a standard ACDM fiat cosmology, with Q m = 0.24, 
» A =0-76, £8 = 0/7 6, h = 0.73 and H = 100 h km s" 1 Mpc" 1 
(ISpergel et alJl2007h . 

2. DATA 
2.1. Object catalogues 

Our galaxies were selected from optical (B W RI) and IR 
imaging (3.6 to 24.0 /mi) in 8.53 deg 2 of the Bootes field 
by the NDWFS, SDWFS and MAGES surveys. Where avail- 
able, we add the NEWFIRM /C-band imaging of 7.80 deg 2 of 
the Bootes field. Any samples selected using the /£ s -band are 
limi ted to this reduced fie ld size. The wide area AGES sur- 
vey dKochanek et al.l2 012) provides redshifts for the brightest 
EROs in the field, while the spectra for the fainter EROs come 
from various small projects using Keck, Gemini and the Kitt 
Peak National Observatory Mayall 4m telescope. 

Source catalogues we re generated by running SExtractor 
dBertin & Arnoutsl [l996) in single-image mode on both the 
/-band and 4.5 /im images. Catalogues were merged af- 
ter accounting for astrometric offsets. We make a pragmatic 
choice to use the 4.5 /im-selected catalogue for source po- 
sitions and photometry, as this catalogue was readily avail- 
able and better represents 3.6 /im sources than /-band. This 
is because objects can be so red as to be detected at 4.5 /im 
but not in the optical. The colors of [3.6] < 18.5 galaxies 
span < [3.6]- [4.5] < 1, so sources with [3.6] = 18.5 will 
have a magnitude of 17.5 < [4.5] < 18.5 and will be com- 
fortably detected in the SDWFS 4.5 /im-selected catalogue 
since it is 80% complete at [4.5] =18.1 and has a 5<r depth of 
[4.5] « 18.5. 

Photometry in all other bands is measured at the 4.5 /im 
detected source positions. For the case of objects detected 
in the 4.5 /im-selected catalogue but with non-detections in 
/-band and /?-band, we assign conservative upper limits of 
/ = 24 and R = 25. These magnitude limits are sufficiently 
bright so that there are no variations in the galaxy distribution 
due to the varying depths of the different pointings. We have 
tested this by inspecting the uniformity of the galaxy distri- 
bution for faint magnitude slices about these assigned limits. 
Our three ERO criteria are shown as applied to the 4.5 /im 
source catalogue (with imaging quality masks applied to re- 
move unusable areas of the Bootes field) in Figure]]] 

We measured aperture photometry using our own analysis 
code, which estimates uncertainties using a Monte Carlo ap- 
proach. The optical and near-IR images were smoothed to a 
common point spread function (PSF) of 1 ."35 (correspond- 
ing to the worst seeing) in the B K RIH K-bands and l."60 in 
/-band, prior to running the photometry code, so that the 
fraction of the flux measured with small aperture photom- 
etry did not vary with position. The fluxes were PSF cor- 
rected and measured using multiple aperture diameters of 
4"- 10" in 1" steps, corresponding to physical distances of 
w 24— 160 /r'kpc at z = 1.0 respectively. We then assign 
the default aperture size for color and magnitude measure- 
ments using a sliding scale based on the 6" aperture /-Band 
apparent magnitude, where fluxes for objects with / > 23.0 
are measured through a 4" aperture and the fluxes for objects 
with / < 20 are measured through 10" apertures, which is 
sufficient to mitigate aperture bias. Astrometric offsets be- 



tween bands were measured and corrected. The residuals 
were fa r smaller than our ape rture size making them neg- 
ligible. Uannuzi & Devi (119991) and iBrown et alf(l2007l) pro- 
vide a dditional details for the NDWFS data and [Ashbv et al. 
(2009) provide details for the SDWFS IRAC photometry. 

The 24 /im fluxes were measured using PSF photometry on 
the MAGES images dJannuzi et al.ll2010b for each source in 
the 4.5 /im-selected catalogue, with uncertainties determined 
using a Monte Carlo method. The NEWFIRM K s catalogues 
are described in Gonzalez et al. (in prep.). We used an ex- 
posure time-weighted map to exclude poor photometry in the 
ZT s -band and we only use photometry from regions where the 
exposure map was greater than 40% of the maximum for the 
NEWFIRM images. Aperture diameters for the NEWFIRM 
K s data were determined as described above. 

To test some of the techniques used in this pa- 
per, we also empl oy 6,285 galax ies observ ed through 
identi cal B H ,RI (iDevetalJ 120071) K s dBundv et al.1 
2006), IRAC (Barmb vetan |2008i) , MIPS (FIDEL; 
bickinson & FIDEL teaml 12007b filters and w ith DEEP2 
DR3 spectra d Yan et all 120091; iNewman et ail 120121 and 
references therein) from the Extended Groth Strip (EGS). 
Photometry for these catalogues was measured and matched 
as described above. 

2.2. Photometric Redshifts 

Photometric redshifts (z p hot) were estimat ed using the ar - 
tificial neural network based ANNz code (iFirth et all 120031) . 
Our spectroscopic tra ining set was largely derived from AGES 
dKochanek et al.l 1201 21) . various small NDWFS p rojects and 
DEEP2 spectroscopy of the EGS (lYan et al.ll2009l) . At z > 0.6 
and / > 20.5, there are relatively few galaxies, so we expanded 
the training set by making fainter copies of bright galaxies 
with modified colors (to preserve color-magnitude relations). 
Objects with colors outside the range of the z < 1 .6 training 
set were not assigned a Zphot, as ANNz has a limited ability to 
extrapolate. Objects with z spe c > 1.6 that are assigned a z P hot 
have systematically underestimated z p hot for the same reason, 
as can be seen in Figure [2] 

Figure [2] demonstrates the accuracy of our z P hot estimation 
for 273 ERO/3.6 galaxies. These galaxies were included in 
the ANNz training set, which allowed it to span the color of 
EROs. We find that 68% of our z P hot estimates for 20 < / < 24 
EROs are within 11% of their corresponding spectroscopic 
redshift (z spec ), or a z /(l+z) ~ 0.06. The z P hot were deter- 
mined using our /-band selected catalogue out to / < 23.5 
and later matched to objects in the 4.5 /im-selected catalogue 
by astrometry. Due to the limited number of z spe c available 
for our ERO s, we cross-ch e cked our z p hot against the tem- 
plate Zphot of Brodw hTet alj (|2006) and find good agreement 
between the two catalogues. 

The fraction of catastrophic Zphot errors, which we de- 
fine as when when A z = \z p hot - z spe c\ > 3<r z , is « 6%. At 
Z w 0.5, stellar spectral features around w 2.3/im produce red 
([3.6]- [4.5]) colors. This can confuse ANNz and is likely re- 
sponsible for the systematic overestimation of z P hot f° r objects 
with 0.4 < Zspec < 0.7 shown in Figure [2] The likely reason 
for systematic z p hot underestimates of objects with z spe c > 1 -6 
is as described earlier. 

We note that our spectra do not uniformly cover the mag- 
nitude and redshift range of our sample. There are very few 
z > 1 .5 galaxy spectra in Bootes and a large number of z < 0.8 
spectra from AGES. Since the very red nature of EROs limits 
the redshift of our sample to z > 0.8, the fraction of z ~ 0.5 
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Figure 1. ERO sample selection from the uncleaned 4.5 ^tm-selected cat- 
alogue. Top to bottom, (R-Ks), (I — K s ) and (/— [3.6]) color for the entire 
[4.5]-selected catalogue in the NDWFS Bootes field. Our ERO selection cri- 
teria (R-K s ) > 5.0, (I -K s ) > 4.0 and (/-[3.6]) > 5.0 are shown by the 
dashed lines. 

galaxies contaminating our sample is probably overestimated 
in Figure |2] We discuss our method for managing these z p hot 
errors, as well as our method for esti mat ing the overall red- 
shift distributions for our samples in 1:14.21 

2.3. Star and Quasar removal 

AGES spectra were used to identify 241 stars and 464 
quasars (QSOs) in our 4.5 /im-selected catalogue. Some of 
the AGES selected stars will have peculiar colors because they 
were targeted as candidate AGN (e.g., X-ray sources). We 
used the (R — I) and (/— [3.6]) color distribution of these stars 
to define a cut to remove these from our ERO samples. This 
star removal cut, of sources with (R—I) > 0.46(7- [3.6])- 
0.26, and with (I - [3.6]) < 5.1, which is shown in Figure [3] 
removes 58,790 objects from the 4.5 ^m-selected catalogue. 
Figure [3] demonstrates that due to the extremely red color of 




Z spec 



Figure 2. Zphot vs - Zspec for 273 ERO/3.6 galaxies in the NDWFS Bootes 
field with both photometric (z p hot) an d spectroscopic redshift (zspec) esti- 
mates. 68% of our Zphot estimates for 20 < / < 24 EROs are within 11% 
(indicated by the dashed lines) of their corresponding spectroscopic redshift 
(Zspec), or cr z /(l +z) 0.06. The results are similar for 20 < / < 22 (trian- 
gles) and 22 < / < 24 (squares), which demonstrate reasonable stability in 
Zphot accuracy between 20 < / < 24. The fraction of catastrophic outliers, 
defined as when A ; = \Zphot —Zspec I > 3fr ; , is 6%. 

EROs, this star cut excludes 100% of stars with a negligible 
impact on the ERO samples. 

There are multi ple IRAC color AGN selection crit eria in 
the literature (e.g., lLacy et alJl2004l IStern et alJl2005h Due 
to the ([3.6] - [4.5]) colors of z > 1 galaxies, the IStern et alJ 
(2006) AGN criterion deliberately removes a large number 
of EROs. Using the colors of AGES identified quasars, we 
demonstrate in Figure [4] the effectiveness of the IRAC color 
selection introduced bv lLacv et al.l (120041) for the ERO sam- 
ples. Figure|4]shows that quasars selected by our ERO criteria 
occupy a r egion of co lor space corresponding to a distinct tail 
wit hin thelLacv et al.l (120041) AGN region. We have modified 
the lLacv et al.l (120041) color criteria to conservatively select 
AGN in this distinct tail seen above the ERO populations in 
Figure [4] It is clear that this selection criteria is not 100% ef- 
fective, but our aim is to select a complete ERO population, 
not a complete sample of AGN. The AGN cut removes 7 1 , 664 
AGN from the 4.5 ^m-selected catalogue. 

After application of both of these stellar and QSO removal 
cuts, our 4.5 /im source catalogue is reduced from 556, 703 to 
426,249 objects. Since the Zf s -band photometry was filtered 
by an exposure time weight map, the 4.5 /im source catalogue 
including good /fj-band photometry is reduced to 366, 070 ob- 
jects. 

3. SAMPLES 

3.1. ERO selection 

From our cleaned 4.5 fj,m source catalogues, the EROrk, 
ERO IK and ERO/3.6 cuts respectively select 16%, 17% and 
27% of all 4.5 ^m detected sources. From 60,066, 63,798 
and 113,098 4.5 /zm-selected ERO RK , ERO !K and ERO /3 . 6 
galaxies, we selected 28,724, 22,451 and 64,370 EROs with 
K s < 19.4, K s < 19.4 and [3.6] < 18.5, respectively. Sam- 



5 




0-[3.6]) 



Figure 3. Control of stellar contamination shown for the uncleaned 4.5 /jm-selected catalogue. (R—I) vs. (/— [3.6]) for All, ERO/jk, ERO/a' and ERO/3.6 objects 
in the NDWFS Bootes field. Green stars are AGES spectroscopically identified stars. We use the location of these stars to place an empirical cut to remove stars 
from our ERO samples. The stellar cut combined with the intrinsically red (/— [3.6]) color of EROs keeps the samples free from stellar contamination. 




([3.6]-[5.8]) 



Figure 4. Control of quasar contamination shown for the uncleaned 4.5 /xm-selected catalogue. A modified form of the Lacy et al. (2004) mid-IR, ([4.5] — [8.0]) 
vs. ([3.6] — [5.8]) AGN selection region is used to minimise quasar contamination within our sample. This is shown for All, EROra:, ERO/k and ERO/3 g objects 
in the NDWFS Bootes field. AGES spectroscopically identified QSOs are shown with red triangles and we use these to p lace a more conservative empirical cut 
at ([4.5] - [8.0]) ^1.8 (solid line), to remove obvious AGN contaminants from our sample. The original Lacy et al. 1 2004) criteria are shown by the dashed lines. 
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pie magnitude limits were set to correspond to the flatten- 
ing in surface density (AO seen in Figure [5] SDWFS is 80% 
complete at [3.6] = 18.2, and this limit is also consistent 
with our 5a /-band and [3 .61 detection limits dJannuzi & Devi 
Il999t lAshbv et all |2009). Throughout this paper we take 
» [3.6] + 1.4, based on the average of K s -[3.6] for our 
EROs. 

The fraction of 4.5 /^m-selected sources meeting the 
ERO/3.6 criterion was found to increase with increasing lim- 
iting magnitude, from 19% at [3.6] < 17.5 to 57% at [3.6] < 
18.5. Similarly, the respective fractions of 4.5 /im-selected 
sources meeting the EROrr- and ERO/a- criteria were also 
found to increase with increasing limiting magnitude from 
11% and 6% at K s < 18.4 to 48% and 35% at K s < 19.4. The 
fractional increase in EROs with increasing apparent magni- 
tude for the EROra- and ERO/A- samples is 29-37%. The 
ERO/3.6 sample has a fractional increase of 38% and selects 
a much larger number of high-redshift galaxies than both 
EROrr- and ERO/a- samples. The EROra- criterion selects a 
greater fraction than the ERO/r- criterion, this is likely due to 
the difference in mean redshift between these two samples. 

As shown in Figure [5] the surface density of the EROrr- 
sample is consistent with that of our previous stu dy of a sub- 
field of Bootes for the same selection criteria (Bro wn et alJ 
12005b . Further more, if we raise our cri teria to (R-K s ) > 5.3 
to match that of Conseli ce et alJ (120081) we find good agree- 
ment in surface density. Likewise, we see good agreement 
between our ERO/y sample and the (I-K) > 4.0 EROs of 
iConselice etaD (120081). To compare with the (i - K) > 3.96 
EROs of iKong et alJ (120091) . we tested the surface density of 
(I—K s ) > 3.7 (which roughly corresponds to (i-K) > 3.96) 
objects and also found good agreement in surface densities. 

In Figure [6] we show the model redshift distributions for 
EROs out to the magnitude limits for each of our samples. The 
mean redshift (I) of the EROra- sample, z = 1.13, is lower than 
that of the ERO/R s ample, with z= 1.19, i n agreement wit h 
previous work (e.g.. IConselice et"al . 2008; Kim et all 1201 ll) . 
We find that the ERO/3.6 sample has a very similar redshift 
distribution to the ERO/R sample, with z = 1.17, although with 
far greater number density. Since all our samples contain at 
least > 9000 objects with photometric redshifts, the statistical 
uncertainties in the mean redshift are negligible (< 1%) and 
we are dominated by systematic uncertainties. 

Figures |3] and |4] illustrate in color space the similarities and 
differences between the ERO/3 6, EROra- and ERO/R samples. 
We see that there is significant overlap between the samples 
selected by the three ERO selection criteria. Quantitatively, 
the bulk (> 86%) of the ERO/r sample is a sub-sample of both 
the ERO/3.6 and EROra- samples. Of the EROra- sample, 78% 
have (/- [3.6]) > 5.0, however only » 55% of the ERO/3.6 
sample have (I-K s ) > 4.0 and (R-K s ) > 5.0. The differences 
are likely due to the ERO/3.6 having a greater surface density 
because at 3.6 we are detecting rest frame NIR emission 
which is less obscured by dust. Thus, both the EROra- and 
ERO/a: samples are largely sub-sets of the ERO/3.6 sample, 
but with different mean redshift. To determine whether there 
is a relationship between ERO clustering strength and lumi- 
nosity, we separate each ERO sample into four [3.6] and K s 
magnitude-limited sub-samples as defined in Tables [5] and [5] 
Again, since all of our samples contain at least > 1000 objects 
with photometric redshifts, the statistical uncertainties in the 
mean redshifts are negligable. 
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Figure 5. Completeness-corrected surface densities (N) for each ERO sam- 
ple. For comparison with the surfac e densities f ound in previous ERO 
studies we show th e distributions from Kone et al. 1 2009) ((i — K) > 3.96), 
IBrownetalj j2005) ((R-K) > 5.0) and Con selice et M JIoOS) ((R-K) > 5.3) 
and ((/ - K) > 4.0). The error bars are included for all points, but in many 
cases are too small to be seen. We use K s [3.6] + 1.4 to shift [3.6] magni- 
tudes onto the K s scale, based on the average of K s — [3.6] for our EROs. 

3.2. ERO Type Separation 

At low redshift the clustering of passive galaxies dif- 
fers f rom that of star-forming galaxies (e.g.. IZehavi et all 
20101). but this may change with redshift (e.g., z > 1.3-2.0; 
Coop er et all 120071 iBrodwin et al.1 120081) . There is no rea- 
son to expect passive and star-forming galaxies to occupy the 
same halos and have identical clustering properties. It is thus 
both interesting and necessary to separate star-forming EROs 
from passive EROs in order to determine if they show differ- 
ing clustering properties. If star-forming EROs have distinctly 
different clustering properties from their passive counterparts, 
it would show that they are two distinct population of galax- 
ies, and it would mean that they are the progenitors of two 
very different populations of local Universe galaxies. It would 
also mean that clustering studies of mixed ERO populations 
cannot be used to accurately link EROs to local galaxy popu- 
lations. 

To independently measure the clustering of passive and 
star-forming EROs, we use the wealth of data in the EGS 
to explore mid-IR c olor separation technique s usin g IRAC 
dWilson et al.1 l2007h and MIPS dStern et al.1 l2006h colors. 
These techniques are physically motivated by polycyclic aro- 
matic hydrocarbon (PAH) emission from warm dust in star- 
forming regions. The ([3.6] - [8.0]) color is able to sepa- 
rate the ERO types as long as the strong PAH emission fea- 
tures lie in 8.0 im \ band, but this ceases to be true for z > 1. 
IStern et al.1 (120061) showed that the difference in ([3.6] - [8.0]) 
color between a star-forming ERO and a passive ERO can 
be very small for archetypal star-forming and passive EROs. 
IKong et al] (120091) also find that very high precision photom- 
etry is required to produce good agreement between ([3.6]- 
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Figure 6. Normalised photometric (dashed histogram) and spectroscopic (solid histogram) redshift distributions as compared to our model distribution (solid 
line) for (from left to right) ERO/3.6 (blue lines) with [3.6] < 18.5, ERO/;k (green lines) with K s < 19.4, and ERO/k (red lines) with K s < 19.4 samples. We find 
that the redder selection criteria (ERO/3 g and ERO/a:) select samples at higher z, as shown by the (z) marker in each panel. While the redshift bins here have 
Az = 0.05 to reduce the noise for the spectroscopic redshifts, Az = 0.02 is used in our analysis. 



[8.0]) color and morphological ERO classifications. We tested 
the ([3.6] -[8.0]) color of EROs in the EGS using DEEP2 
spectroscopic indicators. Although we do observe a tendency 
for passive EROs to have a bluer ([3.6] -[8.0]) color than star- 
forming ones, the two populations are not well separated and 
th is technique is not suitable for this work. 

IStern et all (120061) showed that the (^-[24]) color has a 
much broader color separation than the ([3.6] - [8.0]) color 
for archetypal star-forming and passive EROs. In Figure [7] 
(K s - [24]) color versus K s magnitude is shown for 0.6 < z < 

1 .5 galaxies in the EGS that have B W RIK S , IRAC, 24 fim and 
DEEP2 detections. We independently identify passive EROs 
using the DEEP2 D4000 spectroscopic indicator, which pro- 
vides a measure of the strength of the 4000A break. D4000 
> 1 .5 is used to indicate an evolved stellar population. To in- 
dependently identify the star-forming EROs, we use a DEEP2 
rest-frame [Oil] equivalent width > 5. OA to identify galax- 
ies with significant active star formation. Only DEEP2 spec- 
troscopic indicators with a ^ 3a detection were used in this 
analysis. 

The FIDEL 24 /im catalogue has a 3cr depth of 30 /iJy, 
so we use this as the upper limit on any source without a 
24 fj,m detection. The Bootes field 24 /im data (MAGES) 
has an RMS of 40 /iJy and thus a 2a depth of 80 fjJy. To 
determine how the noisier and shallower MAGES data will 
effect population separation criteria, in Figure [7] we have ap- 
plied additional Gaussian scatter to the FIDEL data to mimic 
the RMS noise level of MAGES. The two populations start 
to merge at K s = 18.4 and [3.6] = 17.5, so the 2a detection 
limit of MAGES prevents us separating the two populations 
at fainter magnitude. 

Since the 3.6 /im imaging in SDWFS is significantly deeper 
than the NEWFIRM K s in Bootes, not every source with a 

3.6 /im detection will have a good K s detection. To achieve 
better separation for the ERO/3.6 sample we tested using the 
([3.6] -[24]) color as a separation method. Figure [8] shows 
that this also works in a similar manner to the (K s — [24]) color. 
In both Figures|7]and[8]the EROs with strong D4000 and those 
with strong [Oil] are well separated. This is further supported 



by the histograms in Figures [7] and [8] which show a bimodal 
distribution in color, but with some overlap, for these two 
spectral indicators. Some galaxies have both strong D4000 
and [Oil]. For these objects we visually inspected the DEEP2 
ID spectra to confirm the color classification and found that it 
accurately distinguished between a passive and star-forming 
galaxies in at least 90% of cases. 

We empirically place a discriminator between passive and 
star-forming EROs where the base of the strong D4000 distri- 
bution (seen in Figures|7]and[8]) and the respective color for an 
object at the 24 /im detection limit meet. Using a 2er detection 
limit for a MAGES source, we can use this separation tech- 
nique out to K s = 18.4 and [3.6] = 17.5, with discriminators at 
(K s - [24]) = 6.0 and ([3.6] - [24]) = 5.2, respectively. We con- 
sider 2a a significant limit for the purposes of this work as we 
are interested in whether the source is emitting at 24 /im rather 
than its specific flux. Since the placement of these discrimi- 
nators is empirical, we varied the discriminators by ±0.2, and 
found no significant change in the separated samples and final 
clustering results. 

We define objects with (K s -[24]) > 6.0 and ([3.6] - [24]) > 
5.2 as star-forming and those with (K s - [24]) ^ 6.0 and 
([3.6] -[24]) 5.2 as passive for K s < 18.4 and [3.6] < 
17.5, respectively. Details for the separated ERO^a-, ERO/^ 
and ERO/3.6 samples are given in Tables [5] and [5] includ- 
ing the sample size, magnitude limit, z and relative fractions 
of the combined ERO samples to the same limiting magni- 
tude. As a check, the archetypal star-forming ERO HR10 
(lGraham&Devl[T996h w ith ffi-[24]) = 7.61 and ([3.6]- 
[24]) = 5.82 (IStern et al.l 120061) is clearly classified as star- 
forming, while th e archetypal passive galaxy LBDS 53W091 
(IStern et al.ll2006l) with (2sT s -[24]) < 5.56 and ([3.6]- [24]) < 
4.13 (assuming a 24 /im flux upper limit of 40 /iJy) is clearly 
classified as passive. 

There is some contamination on both sides of our color 
criteria and the level of contamination increases for fainter 
sources in both Figures|7]and[8] One possible physical reason 
for this is that PAH emission at observed 24 /im varies with 
redshift. There is also rest-frame 10 /im silicate absorption 
which would be observed at 24 /im from a z ~ 1.3 ERO. Con- 
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Figure 7. Galaxy type separation using (K s — [24]) colors as a function of 
K s for all 0.6 < z < 1.5 galaxies in the EGS. Objects with / 24 < 80 /tjy 
are indicated by arrows. Sources are marked by whether they have D4000> 
1.5 indicating an old stellar population (squares), [OII]> 5.0 indicating star 
formation (stars) or both (circles). We empirically separate the populations 
at (K s — [24]) = 6.0. The FIDEL fluxes have been degraded in precision to 
match the noise and detection limit of the Bootes data. The histogram in 
the side panel shows the color distribution of D4000> 1.5 sources (red line) 
and [OII]> 5.0 sources (blue line). The lower panel shows a(K s — [24]) as a 
function of apparent magnitude based on a 24 /xm RMS of 40 fiiy. 

tamination of this nature would therefore increase for fainter 
EROs since these are more likely to have higher redshifts. 

Alternatively, EROs are probably not simply archetypal 
passive and star-forming galaxies. If we assume that the SED 
of a passive galaxy with little dust roughly follows a black- 
body at long wavelengths, we can approximate the expected 
24 /im flux by the Rayleigh-Jeans tail. For LBDS 53W091 
with a 3.6 /im flux of 35 /iJy, we expect an approximate 
24 /im flux of 0.02 /iJy and so should have no 24 fim de- 
tections for truly passive galaxies. However, there are objects 
with both a 24 /im detection and strong D4000 in both Fig- 
ures |7]and|8] The (K s - [24]) and ([3.6]- [24]) colors are ef- 
fectively sampling the specific star formation rate (SSFR) of 
these galaxies. Therefore, these mixed cases are likely passive 
galaxies with some star formation and/or an AGN. To cross- 
check that this is the case, we further examine our cut against 
purely 24 /im non-detected and detected samples, indicating 
star formation rate (SFR) only (discussed in 34.31 1, 

We show the application of these separation criteria to the 
Bootes field EROs/c and ERO/3.6 samples in Figuresl9land[T0l 
respectively. The star-forming EROs form a locus that broad- 
ens toward fainter magnitudes. The broadening is largely an 
effect of the growth in the (K s - [24]) and ([3.6] - [24]) uncer- 
tainties with magnitude, which we show in the lower panel 
of each Figure. Despite the broadening, it is clear that the 
(A^-[24]) and ([3.6] -[24]) colors are finding distinct pop- 
ulations. This is clearly evident in the bimodal distribution 
in ([3.6]- [24]) color shown by the histogram in Figure [TOl 



Figure 8. Galaxy type separation using ([3.6] - [24]) colors as a function of 
[3.6] for all 0.6 < z < 1.5 galaxies in the EGS. Objects with / 24 < 80 tdy 
are indicated by arrows. Sources are marked by whether they have D4000> 
1.5 indicating an old stellar population (squares), [OII]> 5.0 indicating star 
formation (stars) or both (circles). We empirically separate the populations 
at ([3.6] - [24]) = 5.2. The FIDEL fluxes have been degraded in precision to 
match the noise and detection limit of the Bootes data. The histogram in the 
side panel shows the color distribution of D4000> 1.5 sources (red line) and 
[OII]> 5.0 sources (blue line). The lower panel shows <r([3.6] — [24]) as a 
function of apparent magnitude based on a 24 fim RMS of 40 /tjy. 

Figure [9] shows weaker evidence for this bi-modality in the 
(K s - [24]) color, which may be due to the reduced fraction of 
star-forming EROs in the EROrk sample (see, Tables [5] and 
IS. 

The relative fractions of star-forming and passive EROs 
varies between the different ERO selection criteria. At K s < 
18.4, the EROra: sample shows a star-forming fraction of 
30%. Similarly, at K s < 18.4, the ERO/a- sample shows a 
star-forming fraction of 32%. In contrast, at [3.6] < 17.5, 
the ERO/3 6 sample has a significantly higher star-forming 
fraction of 41%. The star-forming fractions found for the 
EROra: and ERO/a: samples are consistent wit h previous stud- 
ies (e .g., the fraction of peculiar galaxies; IConselice et alj 
2008, (R-K s ) > 5.3). The relative fraction of star-forming 
E ROs in the ERO/3 6 sample is consistent with that found 
bvlKong etalJ (120091) using (i-K) > 3.96. iMoustakas eTaT] 
(120041) used a (R-K s ) > 5.0 ERO cut and found a passive 
ERO fraction of only 33-44%, but this study went to a much 
fainter K s < 22 flux limit. We observe in all three ERO sam- 
ples an increase of star-forming ERO fraction with increas- 
ing apparent magnitude (see Tables [5] and |5j and Figures [7] 
[8] and ITOb and this is consistent with the fraction of star- 
forming EROs i ncreasing with decreas i ng luminosity and stel- 
lar mass (e.g., [ Moustaka set alj|2004t : IConselice et alj|2008t 
iKong et alJl2009tlKim et all201 11) . 

4. ANALYSIS 

In this section we analyze the angular and spatial cluster- 
ing of our ERO samples. Using these methods, we place the 
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Figure 9. Separating star-forming and passive EROs using the (K s — [24]) 
color in the EROrk sample. We divide the populations at the dashed line 
(K s — [24]) = 6.0. Upper limits fuu m are marked with arrows. The histogram 
in the side panel shows the relative distribution in color space for Bootes 
objects down to K s < 18.4. The lower panel shows a(K s — [24]) as a function 
of apparent magnitude based on a 24 fim RMS of 40 fi]y. 



EROs in an evolutionar y seq uence and make links to local 
galaxy populations. In § 14. 11 we study the angular clustering 
of each ERO sample as a function of limiting magnitude and 
ERO color and selection criterion. The three ERO samples 
studied, EROrk, ERO/k and ERO/3.6, allow us to explore this 
to some degree across different reds hift intervals. We build 
model redshift distributions in $4.21 and discuss the effects 
and limitations of the photometric redshifts used to generate 
them. We combine the angular clustering and the redshift dis- 
tribution models to estimate the spatial correlation function 
in 34.31 We also analyze the spatial correlation function as a 
function of limiting magnitude and ERO color and selection 
criterion. We use the clustering of EROs and dark matter halos 
to conn ect EROs to other galaxy populations across redshift 
in 34.41 Finally we e stimate the host dark matter halo masses 
of our EROs in j 



4.1. The angular correlation function 

We measured the angular corr elation function for each mag- 
nitude limited sample using the lLandv & Szalavl (fl993l) esti- 
mator, 

A „ DD-2DR + RR 

m= — Tr — , a) 

where, DD, DR and RR are the number of Data-Data, Data- 
Random and Random-Random galaxy pairs at angular sep- 
aration 6± 88/2 respectively. Pairs were counted in loga- 
rithmically spaced angular bins with Alogio 9 = 0.25. Since 
any survey area represents a subset of the Universe, estima- 
tors such as Equation© are subject to an integral constraint 



Figure 10. Separating star-forming and passive EROs using the ([3.6] - [24]) 
color in the ERO/3 5 sample. We divide the populations at the dashed line 
([3.6] - [24]) = 5.2. Upper limits / 24 , i m are marked with arrows. The his- 
togram in the side panel shows the relative distribution in color space for 
Bootes objects down to [3.6] < 17.5. The lower panel shows ct([3.6] - [24]) 
as a function of apparent magnitude based on a 24 /im RMS of 40 /jjy. 



(IGroth & Peeblesl[l977l) . 



w{9 n )dn l dn 2 ~o 



(2) 



where 612 is the angular separation between solid angle el- 
ements dQ\ and dQ.2- Since Equation(0} is normalized by 
itself, the integral constraint leads to a systematic underesti- 
mate of w{6) which is corrected by the addition of the variance 
term, 

1 (3) 



This corrects for the contribution to clustering from the vari- 
ance in number counts over the finite survey area il. We esti- 
mate the uncertainty in each angular correlation function bin 
using a modification to the G aussian approximati on to the co- 
variance matrix discussed in iBrown et al.l (120081) . Our mea- 
sured w(l') values are well below unity, so we consider our 
uncertainty estimates valid to first order. 

The angular correlation functions for the EROrk K s < 18.4, 
ERO ik ^ < 18.4 and ERO /3 6 [3.6] < 17.0 ERO samples are 
presented in Figure [TT] Each of these ERO samples shows 
similar clustering statistics, which is an indication that similar 
galaxy populations are being included by each criteria. We 
then fit the angular correlation function for each sample us- 
ing x 2 statistics including the full covariance matrix for the 
uncertainties and the functional form, 

w(6)=w(l')(^\ \ (4) 
Many studies use a fixed value for 7 (e.g., iFirth et all 
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200llRoche et allllOOl iBrown et al.ll200H iKong et al.ll2009i: 
Kim et alj|201 lh because they have insufficient data to con- 
strain 7. With our larger sample sizes (N > 2000), we are 
able to leave the power law 7 as a free parameter. To bet- 
ter compare with the prior literature and for our smaller sam- 
ples with 1000 < N < 2000, we also fit the correlation func- 
tion w ith a fixed 7 = 1.80 (e.g. jKong et al.l l2009; Ki m et alj 
2011). It should be noted that due to our field size, samples 
with N < 500 do not produce reliable correlation functions 
and have been removed from our final results. 

In many cases the reduced \ 2 for our fits &e distinctly 
greater than unity. We consider this an indication that a sim- 
ple power law is not an adequate description of galaxy angu- 
lar correlation functions. Some studies have tried to compen - 
sate for this with a double power law (e.g.. lKim et alj|201lh . 
However, this is again an empirical approximation and has 
no physical basis. For a robust analysis of galaxy correlation 
we must consider a more complex model of galaxy clustering 
such as galaxy-dark matter halo models (e.g., the Halo Occu- 
pation Distribution; Seliak 2000), which we will address in a 
future paper. 

The angular clustering of EROs is strongly dependent on 
limiting magnitude. As shown in Figure [12] the brightest 
sample for each ERO criterion has a tv(l') clustering ampli- 
tude that is a factor of w 3 stronger than the faintest sample 
(as before, we take K s « [3.6]+ 1.4, based on the average of 
K s — [3.6] for our EROs). A complete summary of our com- 
bined ERO angular clustering analysis results is presented in 
Tables [5] and [5] For ERO samples sub-divided by type, we 
observe stronger angular clustering for passive EROs in all 
three ERO samples. In all cases but the EROrk K s < 17.9 
sample, the vv(l') clustering amplitude of the star-forming and 
passive ERO subsamples straddles that of the combined ERO 
sample. Figure [13] shows angular correlation functions for 
ERO/3.6, [3-6] < 17.5 combined and type-separated samples. 
There is no clear distinction between the passive, star-forming 
and mixed EROra: and ERO/^ samples, but the trend is the 
same as for the ERO/3.6 sample. Separated sample results are 
presented in Tables [5] and[5] 

Freely fitted values for 7 were found to range between 
7 = 1.58- 1.98 for the combined ERO samples with a typ- 
ical uncertainty of 07 = 0.05-0.10. The passive EROs con- 
sistently yield a steeper power law index than the star-forming 
ERO populations (excepting the ERO^- K s < 17.9 sample), as 
can clearly be seen in FigureQ~3] If we take the slope of 7 as 
an indicator of environment, where on the scale of « 1 Mpc a 
steep slope implies a more clustered distribution and a shallow 
7 a smoother distribution, this implies that passive EROs typ- 
ically occupy denser environments than star-forming EROs. 

4.2. ERO redshift distribution 

To transform the angular correlation function to a real-space 
correlation function we need to know the distribution of EROs 
in redshift. The redshift distribution (hereafter dN/dz) for 
each ERO sample was determined by number counts in bins 
of Azphot = 0.02 over < z p hot < 2. The ANNz code as- 
signs galaxies with the same observed photometry to a par- 
ticular z P hot< while the true z S pec for these galaxies will be 
scattered above and below the z P hot and m some cases catas- 
trophically under/over-estimated (see Figure [2). In our par- 
ticular case, since the ERO cut selects a particular redshift 
distribution analogous to a redshift bin, we expect the true 
redshift distribution to be somewhat broader than given by 
the dN/dz (Zphot)- This effect is clearly observed in Figure[6] 



We therefore need to model the true dN/dz (z) based on the 

dN/dz (Zphot)- 

We build our dN/dz (z) models in a two stage process. 
Firstly, we replace the z p hot of each galaxy with a Gaussian 
probability distribution accounting for 90% of each galaxy's 
distribution in redshift, where = Zphot and c is a function of 
/-band magnitude, {Bw~R) color and Zphot- This accounts 
for the accuracy of the photometric redshift itself. Secondly, 
we account for catastrophic Zphot errors with a second Gaus- 
sian probability distribution making up the remaining 10% 
of each galaxy's distribution in redshift, where /i = z P hot and 
a = 0.5. This corresponds to the ps 6% of galaxies having 
A- > 3<7 z z/(l +z), as observed i n the data. The outliers need 
not be precisely modelled, as the Limber] d 19541) equation (see 
34.3b depends on (dN/dz) 2 - The choice of 10% was made 
to better match the z spe c distribution. Our resultant dN/dz (z) 
models are reasonable representations of the true redshift dis- 
tribution, as shown by the examples in Figure [6] 

The z for each of our ERO samples is giv en in Tables [5] 
[5] [5] and [5] and for the reasons discussed in 32.21 we do not 
quote uncertainties on these values. As demonstrated by Fig- 
ure[6] the EROrk sample has lower z than both the ERO/^ and 
ERO/3.6 samples, while the ERO/3.6 and ERO«: have similar 
Z. Fainter magnitude limited samples have higher z, which is 
expected as a fainter limiting magnitude will include a larger 
number of more distant galaxies. The passive EROs and star- 
forming EROs in the type-separated samples have very simi- 
lar z to each other. 

Reliable z p hot are only assigned to objects with / < 23.5, 
therefore not all of the EROs in our sample are assigned a 
Zphot- We account for this in two ways. For EROs with no 
Z p hot> we assume that they have a dN/dz (z) similar to those 
with Zphot < and simply scale up dN/dz (z) to the full sample 
size. We do this for completeness, and it has li ttle impact 
on the clusteri ng re sults because the iLimberl (119541) equation 
(discussed in 34.31 ) only depends on the shape of dN/dz (z). 
We also provide in Tables[5]and[5]a complete clustering anal- 
ysis including only objects with z p i 10 r information for compar- 
ison. In all but our faintest magnitude limited samples, > 50% 
of EROs in each sample have z p hot information and the mod- 
elled dN/dz is considered a valid representation of the entire 
sample. 

It should be noted that when using smaller sample sizes 
but with an identical dN/dz (z), the clustering length of the 
Zphot only samples (shown in Tables [5] and [5]), is observed to 
be systematically larger than that of samples including EROs 
without Zphot information (shown in Tables [5] and [5]). For the 
majority of samples, the increase in ro is < 20%, with the ex- 
ception of samples that have a large number of objects without 
Zphot information. 

We speculate that this is because objects without Zphot in- 
formation are likely faint and therefore more likely to have a 
higher average redshift. When included in the angular clus- 
tering measurement, these high-z objects increase the number 
of uncorrelated pairs and will reduce the angular clustering 
amplitude. This would be accounted for if the redshift distri- 
bution was broader to account for the higher redshift of these 
faint objects. We note that samples with and without com- 
plete Zphot information show the same general trends. Since 
we get better clustering statistics with greater sample size, we 
use our samples including EROs without Zphot f° r discussion 
and the samples using only objects with a measured Zphot are 
included as a cross-check. 




Figure 11. Angular correlation functions w(9). From left to right we show w(9) for the ERO/3.6 [3.6] < 17.0, EROrk K s < 18.4 and ERO/k K s < 18.4 samples. 
We also show the reduced x 2 an d fitted power law 7 for each fit. These three samples are selected using different ERO criteria but show very similar clustering 
statistics. This suggests similar populations are being selected by all three ERO criteria. 
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Figure 12. The angular correlation strength of our ERO samples as a func- 
tion of apparent magnitude and for a fixed slope of 7 = 1.80. For clarity, mag- 
nitude offsets of +0.05 and —0.05 were applied to the ERO/3 6 and ERO;^ 
data, respectively. Angular clustering is seen to be strongly dependent on 
sample limiting magnitude. Here we have used K s ss [3.6] + 1.4, based on the 
average of K s — [3.6] for our EROs, to place the [3.6] samples on the figure. 

4.3. Spatial clustering 

The spatial cor relation function £(r) is related to w{9) by 
the iLimbeil (1 1 9541) equation, 



tv(0): 
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The characteristic clustering length ro gives a measure of the 
clustering strength, and is related to £(r) by its functional form 



(6) 
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where ro is a measure of the spatial separation at which we ex- 
pect to find an excess pair of galaxies above a random galaxy 
distribution, and 7 is as previously measured from the angu- 
lar correlation function. The results of our spatial clustering 
analysis are given in Tables [5] and [5] for the combined ERO 
samples and in Tables [5] and [5] for the separated ERO sam- 
ples. The uncertainty in ro was determined from the spread 
in ro values of the angular correlation function models with 

x 2 <(xL+i)- 

Using our magnitude-limited samples we explore the spa- 
tial clustering of EROs as a function of limiting luminosity. 
For clarity and ease of comparison to previous studies we 
will only discuss results using fixed 7 = 1.80. In Figure [141 
it is clear that ERO clustering strength is a strong function 
of apparent magnitude, since EROs occupy a particular red- 
shift range, magnitude is a rough proxy for stellar mass. The 
EROra: and ERO«: samples show stronger clustering than the 
ERO/3.6 samples at similar magnitude limits (as before, we 
take K s sa [3.6] + 1 .4, based on the average of K s - [3.6] for 
our EROs). This could be due to the ERO/3.6 sample having a 
larger star-forming galaxy fraction. 

For the ERO«a: and ERO/a: samples at a fixed magnitude 
limit, the clustering strength depends o n the z of each s ample, 
in agreement with previous work (e.g.. lKim et al.ll20Tlh . This 
difference is easily accounted for because the difference in z 
between these two samples makes the objects in the ERO/a: 
sample on average s» 6% brighter for a given limiting magni- 
tude than the EROra: sample. With this in mind we speculate 
that there is very little evolution in ro for roughly equivalent 
mass EROs between z = 1.2 and z= 1.0. The difference be- 
tween the ERO/3.6 sample and the EROra: and ERO/a: sam- 
ples, we attribute (as above) to the difference in relative star- 
forming to passive ERO fractions. 

The brightest ERO /3 . 6 , EROra- and ERO/a- EROs have clus- 
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Figure 13. Angular correlation functions w(9) for (from left to right) the combined, passive and star-forming EROs brighter than [3.6] < 17.5 selected with 
(/— [3.6]) > 5.0. We show the reduced \~ an d fitted power law 7 for each respective fit. Passive EROs show significantly stronger angular clustering strength 
than those with star-formation, especially on small angular scales. 



tering lengths of r = 10.14 ±0.76 /T 1 Mpc, r = 10.73 ± 
0.88 h~ l Mpc and r = 9.88 ± 1.54 /r 1 Mpc, respectively, 
comparable to tha t of « 4L* local ellipticals (SDSS; e.g., 
IZehavi et al.1l2010t) . The faintest of ERO RK and ERO, K sam- 
ples are still strongly clustered, with ro = 6.49 ± 0.34 h Mpc 
and ro = 7. 10 ±0.36 h~ l Mpc, respectively. This is com- 
parable to that of 2L* local ellipticals. With a significantly 
higher star-forming fraction, the faintest ERO/36 sample 
shows ro = 4. 90 ±0.20 h~ l Mpc, similar to w L* local galaxies 
(SDSS; e.g.. IZehavi et al.ll2010t) . It is clear that a larger star- 
forming fraction leads to reduced clustering strength. Since 
the star-forming fraction increases with apparent magnitude, 
this likely contributes to the observed decline of the combined 
ERO sample clustering length with apparent magnitude seen 
in Figure [14] 

The strong clustering of the faint ERO«a- and ERO/^ sam- 
ples with respect to the faint ERO/3.6 samples suggests that 
there is a dependence of clustering strength on ERO type. 
The clustering results for our separated ERO^- and EROik 
K s < 17.9 samples are largely inconclusive due to statistical 
uncertainties. However, all the type-separated ERO samples 
show a difference between the clustering of passive and star- 
forming EROs. In the ERO/3 6 samples, which have the small- 
est uncertainties, there is a stark difference in the clustering 
between the two ERO types, with the clustering length ro of 
the passive EROs being a factor of w 1 .4 greater than that for 
the star-forming EROs at all 3 magnitude limits. It is clear that 
star-forming EROs occupy very different environments from 
passive EROs with co mpara ble apparent magnitudes, and this 
is discussed further in 34.41 

As a cross-ch eck for the reliability of our ERO separa- 
tion technique ( 33.21 ). we tested in the ERO/3 6 sample out 
to [3.6] < 17.5, the clustering of the 9,449 24 fim non- 
detections. For a fixed 7 = 1.80 we estimate ro = 8.46 ± 
0.48 /r'Mpc, and this agrees well with the clustering length 
of the corresponding passive sample. We then tested the clus- 
tering of the 12,344 24 /im detections for a fixed 7 = 1.80 
and estimate ro = 6.22 ± 0.40 /r'Mpc, which is similarly in 
agreement with the clustering length of the corresponding 



star-forming sample. Both results are slightly higher than for 
the color separated EROs. We take this as evidence that color 
classified passive EROs with 24 /im detections increase the 
measured clustering length of the sample with 24 /im detec- 
tions, while star-forming EROs are absent in the sample with 
no 24 /im detections and hence there is less likelihood of con- 
tamination. This conclusion is supported by the majority of 
color classified passive EROs with 24 /mi detections having 
strong D4000, as seen in Figures [7] and [8] This confirms ERO 
type definitions based on 24 /im colors as a robust method for 
separating star-forming and passive EROs. 

4.4. Linking EROs throughout z via spatial clustering 

The clustering of passive and star-forming EROs with com- 
parable apparent magnitude show that the two ERO sub- 
populations occupy very different environments. This sug- 
gests that these sub-populations of EROs are the progeni- 
tors of different local Universe galaxies. In order to make 
links between EROs and galaxy populations at other epochs, 
we present comparisons of clustering for passive and star- 
forming [3.6] < 16.5 EROs with samples of some of the 
largest galaxy clustering studies over < z < 2.5. 

For ERO/3.6 [3.6] < 16.5 passive EROs we use the very 
brightest samples from each available comparison study and 
place our results in this context in Figure [15] We have also 
determined the expected ro(z) for dark matter halo samples 
and a range of m inim um halo masses based on the methodol- 
ogy outlined in 34.51 The clustering at fixed minimum halo 
mass declines with redshift for massive halos (e.g. Mhaio > 
10 13 Mq). Halo masses however, grow with time and we 
would expect ro for a single, massive, growing halo to be 
slightly increasing with decreasing redshift. This is also the 
case for the massive galaxies at their cores. Figure [TBI shows 
that reasonably bright passive EROs are the likely prog enitors 
of > 4L* local Universe ellipticals (IZehavi et alil2010i) . This 
suggests that bright EROs are the likely desce ndants of lumi- 
nous dust-obscured galaxies (DOGs) at z ~ 2 (B rodwin et alj 
2008) and the likely progenito rs of some of the m ost massive 
galaxies in the local Universe (Ze havi et all 2010). 
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Figure 14. The real space spatial clustering strength of EROs as a function of 
apparent magnitude with fixed 7 = 1.80 for the three ERO samples. For clar- 
ity, magnitude offsets of +0.05 and -0.05 have been applied to the ERO/3 g 
and ERO/k data, respectively. The spatial clustering is strongly dependent on 
sample limiting magnitude. Here we assume K s ss [3.6] + 1.4, based on the 
average of K s — [3.6] color of our EROs. 
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Figure 15. Spatial clustering as a function of redshift for the passive ERO73.fi 
[3.6] < 16.5 sample (red circle; with fixed 7 for ease of comparison). The 
IZehavi et al] 120101) . Mf" = -22.0 sample is indicated by a triangle, the 
IWake et al] J201 If ) , stellar mass of = 10 11 Mq sample with a diamond and 
the most luminous Brodwin et al. 1 2008) dust obscured galaxy sample by the 
blue square. The curves show the evolution of rg(z ) for fixed minimum halo 
mass, based on the methodology discussed in 84.51 



Figure 16. Spatial clustering as a function of redshift for the star-forming 
EROy3 5 [3.6]<17.5 sample (red circle; with fixed 7 for ease of comparison). 
Here we compare to the brightest absolute magnitude limited samples of pas- 
sive red galaxies from Brown et al. 12008), for eac h redshift range (g reen 
stars). The most massive, stellar mass samples from Wake et al. 12011), for 
each redshift range (diamonds), a nd the least luminous dust obscured galaxy 
sample from Brodwin et al. (2008) (blue square). The curves show the evo- 
lution of r p(z) for fixed minimum halo mass, based on the methodology dis- 
cussed in j )4.5l 

In the case of our ERO/3.6 [3.6] < 16.5 star-forming EROs, 
their clustering strength of ro = 8.89 ± 1.73 h ~ l Mpc, asso- 
ciates them with > 2L* local red galaxies (e.g., Zeh avi et alj 
2010). Furthermore, when placed into context with compar- 
ison studies (Figure [TBI), the expected evolution in clustering 
for the ERO/3.6 [3.6] < 16.5 star-forming EROs makes these 
likely progenitors of z < 0.8 passive galaxies. Furthermore, 
the clustering length of [3.6] < 16.5 star-forming EROs is 
comparable to that of passive [3.6] < 17.5 EROs. This in- 
dicates the possibility of a duty cycle, whereby star-forming 
EROs turn into passive EROs with significantly fainter ap- 
parent magnitudes, which builds a consistent picture. This 
is evidence that > L* local red galaxies are building up stellar 
mass at z ~ 1.2. The 24 /im flux limit of MAGES is consistent 
with tha t of Luminous Infra-R ed galaxy (LIRG) emission at 
Z « 1.2 dLe Floc'h et al.l l2005). whic h sets the minimum star 
formation rate of at rj 10 M Q year" 1 dNoeske et al.ll2007l) for 
star-forming EROs with [3.6] < 17.5. These galaxies may be 
undergoing merger-induced star formation or are early forms 
of present-day ellipticals with continuing/residual star forma- 
tion. 

The above associations provide a strong indication that the 
ERO population is an earlier form of the elliptical population 
we observe in the local Universe. In general, the ro values 
for star-forming and passive EROs straddle that of the overall 
ERO population of the same apparent magnitude, with passive 
EROs being more strongly clustered than star-forming EROs 
for a particular magnitude limit. 

4.5. Host halo mass 

As galaxies are biased tracers of dark matter, it is possi- 
ble to get an estimate of their host halo mass through a direct 
comparison of their clustering with that of dark matter ha- 
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los. We determine the clustering of dark matter halos from 
the non-linear dark matter power spectrum. The power spec- 
trum of dark mat t er is g enerated using the transfer function of 
l Eisenstein & Hul (119981) . the halo mass function and bias of 
Isheth & Tormenl (1 19991) . and the halo model of lSeliakl (120001) 
We a ssume NFW dark matter halo p rofiles (iNavarro et al] 
1996) with the concentration scalings of Bull ock et al.l ( 2001 ). 
We adopt a simple halo occupation scheme assuming one 
galaxy per halo. We assume that the EROs represent the spa- 
tial distribution of their host halos on the scales of «Mpc, and 
determine their host halo mass by matching the ro of the EROs 
with that of mass limited dark matter halo samples. 

This scheme provides crude estimates of the host halo mass, 
as it neglects the possibility of ERO satellite galaxies, and 
the correlation function of galaxies is not well described by 
a power law. A more accurate treatment may be achieved 
by halo occupation distr i bution modeli n g (e .g., Selia kil2000b 
iPeacock & Smith! 120001: IZhengl l2004t IZehavi et all 12004 
which we will address in a future paper. We show the cor- 
responding halo mass estimates for combined ERO samples 
in Tables [5] and [5] and the separated ERO samples in Tables 
and |5] We find our brightest EROs are associated with dark 
matter halos with masses of the order 10 13 M Q . This host halo 
mass agrees with that predicted from the halo occupation dis- 
tribution by Gonzale z-Perez et"afl (1201 ll) . The faintest EROs 
are associated with dark matter halos of m 10 I2 M Q . These 
halo masses agree with the ERO-local galaxy population as- 
sociations made using the spatial clustering measurements in 
the previous section. 

5. DISCUSSION & SUMMARY 

We present the most statistically robust measurement of 
ERO clustering to date. The effects of cosmic variance on 
our results are modest considering that the integral constraint 
increases w(9) by only « 3% due to our large survey area. Our 
redshift distribution was derived from photometric redshifts, 
which are constrained with spectroscopy, convolved with their 
photometric uncertainties, and account for catastrophic z P hot 
outliers, making them reliable representations of dN/dz (z). 
We find clear trends of decreasing clustering strength with 
apparent magnitude, which roughly translates to stellar mass 
since the redshifts are similar. Previous studies have hinted at 
such a relationship, but t he samples were too small to draw 
solid conclusions (e.g.. iBrown et alj|2005l) . Our EROrk re- 
sults agree with our previous m easurement o f ERO clustering 
in a sub-field of the NDWFS (Br own et al.ll2005l) . We find 
that our measured ERO^ and ERO/a: clustering lengths are 
typically l ower than those of t he previous l argest area ERO 
studies bv lKim et all d2011l) and lKong etail (120091) . but agree 
within 2a. It shou ld be noted that both lKim et aiTd201 lb and 
iKong et ail (120091) have likely underestimated their uncertain- 
ties, as they are based on Poisson uncertainties and assume a 
diagonal covariance matrix, so the true level of disagreement 
is smaller. 

We tested whether the difference in ERO clustering results 
can be explained by the different criteria used to define ERO 
samples. We have clearly shown that the various ERO crite- 
ria select similar but not identical galaxy populations. Differ- 
ences in number density, z and star-forming to passive ERO 
fraction are likely the main causes for differing estimates of ro 
between selection criteria. The results of ERO clustering stud- 
ies using different selection criteria should be directly compa- 
rable so long as these factors are considered. For example, 
ERO ik K s < 18.4 EROs possess slightly stronger clustering 



properties to EROrk K s < 18.4 EROs, because the difference 
in z between the samples means that the ERO/a- EROs are 
systematically brighter than the ERO«a: EROs. In the case of 
ERO/3.6 we found that this criterion selected most of the same 
galaxies as the EROa'a' and ERO/a: samples, but also a greater 
surface density of objects for a comparable magnitude. The 
additional objects included a larger fraction of star-forming 
EROs than for the EROrk and ERO/a- samples, and this sam- 
ple then showed systematically lower clustering for similar 
magnitude limits. This means that when you are trying to link 
EROs to local galaxy populations by their clustering, you can- 
not obtain meaningful results unless you separate the passive 
and star-forming EROs into separate samples. 

Our separation criteria, based on the work of Ster n et ail 
(2006), shows strong separation in both (K s — [24]) and ([3.6]- 
[24]) colors between EROs with star formation and those 
without. The success of this technique is its ability to pick 
out PAH emission from hot dust at high-re dshift and oyer the 
redshift range re quired for ERO samples dStern et al J 120061 : 
iTrever et al.ll201of) . The ([3.6] -[8.0]) color is less effective 
because 8.0 /im is no longer sampling strong PAH features at 
Z > 1. Likewise the {J-K) color is not useful for most EROs, 
because it does not bracket or straddle the 4000A break until 
Z > 1.4. Deeper 24 data would allow this technique to be 
used for fainter samples than analysed in this work. 

At fixed apparent magnitude, we find that the two ERO 
sub-populations occupy different environments, with passive 
EROs being more strongly clustered than those with star for- 
mation. The clustering strength for each type-separated ERO 
sample is also a function of apparent magnitude. In many 
cases, passive EROs have similar clustering properties to star- 
forming EROs of brighter apparent magnitude. This is consis- 
tent with a duty cycle, whereby star-forming EROs turn into 
passive EROs with significantly fainter apparent magnitudes. 

The SSFR seems to be lower in dense environments at low 
redshift (e.g. IZehavi etall 1201 Oh . but s ome studies sugges t 
that this tren d changes at z > 1.3 (e.g.. ICooper et al.ll2007l) 
to z w 2 (e.g jBrodwin et al.ll2008tlMaghocchetti et al.112 008). 
While we observe that EROs without star formation typically 
occupy denser environments than star-forming EROs, we see 
evidence for star formation in much denser environments than 
observed at z < 1. This suggests that a transition between 
over-densities being dominated by passive galaxies to being 
dominated by star-forming galaxies coul d well be occurring at 
z > 1 .3, in agreement with the findings of lCooper et al.l (2007) 
and lBrodwin et al.l (120081) . 

Denser environments imply higher mass dark matter ha- 
los, and star-forming EROs seem to have host halo masses 
of up to M/ W i « 10 13 M Q , an order of magnitude greater 
than t he host halo mass of a typical star-forming galaxy at 
z = {Zheng et al. 2007). The shifting of the peak in star for- 
mation to higher mass halos at z > 1 is in agreement with 
the results of previ ous hi gh-redshift c lustering studies (e.g., 
Magli occhettietaLl 12001 [Wake et al.ll20TTh . This could be 
evidence for "downsizing", where the host halo mass and 
stellar mass of star-forming galaxies decreases with decreas- 
ing redshift. Since a fixed halo mass at high redshift is 
more bias ed with respect to the underlying matter distribu- 
tion (e.g.. ISheth & Tormenl [1991 ISeliak[l2000l) and therefore 
has a greater merger and mass accretion rate than at z < 1, this 
result is not unexpected. 

The most significant result of our study is that the pas- 
sive and star-forming ERO sub-populations are two distinct 
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galaxy populations. Separating the ERO population into its 
sub-populations has allowed us to place EROs into an evolu- 
tionary sequence. Bright passive EROs are the likely progen- 
itors of local > AL* red galaxies, indicating that massive red 
galaxies have had their star formation truncated prior z ~ 1 .2. 
We were not able to obtain a clustering measurement for the 
progenitors of central cluster galaxies as such structures are 
rare, and as a consequence there are too few of these galax- 
ies to measure their clustering. However, our brightest ERO 
samples certainly hint that the very brightest EROs are the 
likely progenitors of the most massive ellipticals in the local 
Universe. 

Interestingly, the clustering of star-forming EROs connects 
them as the likely progenitors of local > L* red galaxies. It 
may seem pec uliar that we detect significant star formation 
(> 10 M /yr; iNoeske et aTll2007l) at z « 1.2 in progenitors 
of local > L* ellipticals. As it is well established that star for- 
mation is expected to have ceased in these objects at zk 1.5 
(e.g., gravitational lensing; fRusin & Koc hanek 2005). How- 
ever, red/elliptical galaxies can be modelled with tau models 
with e-folding times of t < lGyr and formation redshifts of 
Z > 2. It is thus possible to have a galaxy where the bulk of 
star formation took place at z > 1.5, but at z= 1.2 the SFR 
is still ss 10 M Q year" 1 while by z = 0.5 the SFR is negligi- 
ble. Therefore, we find that local red galaxies, down to L*, in 
some cases were still building stellar mass via star formation 
at this epoch. 
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Table 1 

Summary of combined ERO clustering results 





ERO cut 


Mag. limit 


ERO type 


N 


z 


7 


w(l') 


X 2 /d.o.f 


r (h 1 Mpc) 


M halo (xlQ 12 M Q ) 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


ALL 


3209 


1.11 


1.91 ±0.08 


0.295 ±0.028 


8.48 


8.77 ±1.65 


6.725 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


ALL 


9537 


1.14 


1.98 ±0.06 


0.182±0.013 


2.29 


6.78 ±0.95 


2.9 + \i 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


ALL 


21800 


1.16 


1.89 ±0.05 


0.126 ±0.008 


2.93 


6. 10 ±0.77 


2.1 + J : » 


I- 


[3.6] > 5.0 


[3.6] < 18.0 


ALL 


40969 


1.17 


1.87 ±0.05 


0.085 ±0.005 


9.02 


5.08 ±0.54 


1 1+ 04 

-0 4 


I- 


[3.6] > 5.0 


[3.6] < 18.5 


ALL 


64370 


1.17 


1.76 ±0.02 


0.071 ±0.002 


22.92 


4.93 ±0.21 




R 


-K s > 5.0 


Ks < 17.9 


ALL 


2715 


1.05 


1.77 ±0.07 


0.312±0.031 


2.95 


11.21 ±2.08 


14.3!« 


R 


-K s > 5.0 


Ks < 18.4 


ALL 


6801 


1.09 


1.80 ±0.06 


0.238 ± 0.020 


3.86 


9.63 ±1.57 


8.9!" 


R 


-K s > 5.0 


K s < 18.9 


ALL 


15013 


1.12 


1.75 ±0.06 


0.169 ±0.014 


3.46 


8.43 ±1.39 


5.9 + l« 


R 


-K s > 5.0 


Ks < 19.4 


ALL 


28724 


1.13 


1.83 ±0.06 


0.1 12 ±0.008 


7.34 


6.26 ±0.93 


2-3^ 


I 


-Ks > 4.0 


Ks < 17.9 


ALL 


1314 


1.08 


1.58±0.10 


0.235 ±0.041 


2.80 


12.68±3.14 


21.0 + 2§o 


I 


-K, > 4.0 


Ks < 18.4 


ALL 


3672 


1.15 


1.78 ±0.08 


0.295 ±0.029 


3.85 


11.19±2.27 


14.2+^° 


I 


-Ks > 4.0 


Ks < 18.9 


ALL 


9829 


1.18 


1.71 ±0.07 


0.223 ± 0.022 


6.92 


9.98 ±2.00 


10.0+™ 


I 


-Ks > 4.0 


Ks < 19.4 


ALL 


22451 


1.19 


1.75 ±0.07 


0.148 ±0.014 


9.83 


7.62 ±1.35 


a 9+2.6 
-9 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


ALL 


3209 


1.11 


1.80 


0.319 ±0.024 


8.21 


10. 14 ±0.76 


10.5+^ 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


ALL 


9537 


1.14 


1.80 


0.213 ±0.012 


3.09 


8.53 ±0.48 


6.1-} i 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


ALL 


21800 


1.16 


1.80 


0.137 ±0.007 


2.73 


6.81 ±0.35 


3-0+J! 


I- 


[3.6] > 5.0 


[3.6] < 18.0 


ALL 


40969 


1.17 


1.80 


0.092 ±0.004 


7.69 


5.53 ±0.24 


1 5+ ' 2 


I- 


[3.6] > 5.0 


[3.6] < 18.5 


ALL 


64370 


1.17 


1.80 


0.073 ±0.003 


13.90 


4.90 ±0.20 


1 0+ 01 

1 -"-o.i 


R 


-K s > 5.0 


K s < 17.9 


ALL 


2715 


1.05 


1.80 


0.305 ± 0.025 


3.02 


10.73 ±0.88 


12.5+H 


R 


-K s > 5.0 


K s < 18.4 


ALL 


6801 


1.09 


1.80 


0.238 ±0.015 


3.86 


9.63 ±0.61 


8 9+ 1 ' 


R 


-K s > 5.0 


Ks < 18.9 


ALL 


15013 


1.12 


1.80 


0.161 ±0.008 


3.94 


7.86 ±0.39 


4 7+0.8 


R 


-K s > 5.0 


K s < 19.4 


ALL 


28724 


1.13 


1.80 


0.1 15 ±0.006 


7.03 


6.49 ±0.34 


2 9+ L5 


I 


-Ks > 4.0 


Ks < 17.9 


ALL 


1314 


1.08 


1.80 


0.225 ±0.035 


3.54 


9.88 ±1.54 


9.732 


I 


-Ks > 4.0 


Ks < 18.4 


ALL 


3672 


1.15 


1.80 


0.291 ±0.021 


3.94 


10.89 ±0.79 


13.1+1;' 


I 


-Ks > 4.0 


Ks < 18.9 


ALL 


9829 


1.18 


1.80 


0.205 ±0.011 


8.72 


8.77 ±0.47 


6 7+ L2 


I 


-Ks > 4.0 


Ks < 19.4 


ALL 


22451 


1.19 


1.80 


0.140 ±0.007 


11.85 


7.10±0.36 


T. 4+0.6 
J -+-0 4 



Note: Clustering results for samples with N < 500 are not shown due to poor quality statistics for our field area. 



Table 2 

Summary of type-separated ERO clustering results. 





ERO 


Mag. 


ERO 


N 


% 


z 


/ 




Y 2 Id o f 


' U 


lvl nalo 




cut 


limit 


tvDe 




total 










(/r 1 Mpc) 


(xlO 12 Ma,) 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


E 


2281 


71 


1.11 


1.92 ±0.09 


0.419 ±0.041 


5.91 


10.01 ± 1.79 


10. l +f ?S 
" -4.7 


1- 


[3.6] > 5.0 


[3.6] < 16.5 


SF 


928 


29 


1.09 


1.82±0.14 


0.232 ±0.045 


7.83 


8.75 ±2.27 


6.6+1? 

-4.1 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


E 


6183 


65 


1.14 


2.00 ± 0.07 


0.265 ±0.020 


4.17 


7.73 ± 1.05 


4 5 + ?i 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


SF 


3354 


35 


1.14 


1.83 ±0.10 


0.139±0.016 


6.53 


7.05 ± 1.45 


3.3+?^ 

-1.7 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


E 


12905 


59 


1.16 


1.95 ±0.06 


0.185 ±0.013 


2.98 


6.84 ±0.88 


3.0 + | f 

— 1.1 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


SF 


8895 


41 


1.16 


1.87 ±0.09 


0.091 ±0.008 


5.52 


5.52 ± 1.02 


1.5+V' 

-0.7 


R 


-K s >5.0 


Ks < 17.9 


E 


2213 


82 


1.04 


1.89 ±0.08 


0.335 ±0.033 


3.02 


10.36 ± 1.87 


-4.9 


R 


-K s > 5.0 


Ks < 17.9 


SF 


502 


18 


1.06 


2.04 ±0.17 


0.401 ±0.075 


4.50 


9.20 ±2.70 


7.7 + ^ 

-5.2 


R 


-Ks > 5.0 


K s < 18.4 


E 


4744 


70 


1.08 


1.85±0.06 


0.260 ± 0.022 


5.56 


9.52 ± 1.50 


8.6+i? 

-3.5 


R 


-K, > 5.0 


Ks < 18.4 


SF 


2057 


30 


1.11 


1.79 ±0.09 


0.262 ±0.031 


1.97 


10.34 ±2.23 


11.1+9.3 

-5.5 


I 


-K, > 4.0 


Ks < 17.9 


E 


1114 


85 


1.07 


1.73±0.12 


0.247 ± 0.042 


2.37 


11.37 ±2.90 


14.9+' 5 , 2 

-8.8 


I 


-K s > 4.0 


K s < 17.9 


SF 


200 


15 


1.15 












I 


-Ks > 4.0 


Ks < 18.4 


E 


2481 


68 


1.13 


1.78 ±0.09 


0.313 ±0.033 


5.14 


11.60 ±2.46 


15.9 + ' 2 „" 

-7.8 


I 


-Ks > 4.0 


Ks < 18.4 


SF 


1191 


32 


1.18 


1.61 ±0.10 


0.332 ±0.051 


3.27 


13.95 ±3.33 


28.1+?|-2 

-IS. 7 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


E 


2281 


71 


1.11 


1.80 


0.456 ±0.035 


6.00 


11.82±0.91 


16.9+" 

-3.7 


1- 


[3.6] > 5.0 


[3.6] < 16.5 


SF 


928 


29 


1.09 


1.80 


0.231 ±0.045 


7.76 


8.89± 1.73 


6.9+H 

-3.4 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


E 


6183 


65 


1.14 


1.80 


0.315±0.019 


4.86 


10.12±0.61 


" -1.8 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


SF 


3354 


35 


1.14 


1.80 


0.141 ±0.016 


6.38 


7.27 ±0.83 


3.7+J-i 

-1.2 


I- 


[3.6] >5.0 


[3.6] < 17.5 


E 


12905 


59 


1.16 


1.80 


0.212±0.011 


3.43 


8.29 ±0.43 


5.6+10 

— u.y 


I- 


[3.6] >5.0 


[3.6] < 17.5 


SF 


8895 


41 


1.16 


1.80 


0.095 ±0.007 


5.39 


5.91 ±0.44 




R 


-K s >5.0 


K s < 17.9 


E 


2213 


82 


1.04 


1.80 


0.357 ±0.030 


3.01 


11.77 ±0.99 


16.(^1 


R 


-K s >5.0 


K s < 17.9 


SF 


502 


18 


1.06 


1.80 


0.400 ±0.080 


4.42 


11.69 ±2.34 


16.3+^ 2 2 ' 


R 


-Ks >5.0 


Ks < 18.4 


E 


4744 


70 


1.08 


1.80 


0.271 ±0.019 


5.41 


10.23 ±0.72 


10.8+ 2 -5 


R 


-K s >5.0 


K s < 18.4 


SF 


2057 


30 


1.11 


1.80 


0.261 ±0.027 


2.00 


10.21 ±1.06 


io.7+3 : 8 


I 


-Ks > 4.0 


Ks < 17.9 


E 


1114 


85 


1.07 


1.80 


0.246 ±0.040 


2.52 


10.56 ±1.72 


11.9+^1 


I 


-Ks > 4.0 


Ks < 17.9 


SF 


200 


15 


1.15 












I 


-Ks > 4.0 


Ks < 18.4 


E 


2481 


68 


1.13 


1.80 


0.309 ±0.027 


5.20 


11.29±0.99 


14 6+ 41 


I 


-Ks > 4.0 


Ks < 18.4 


SF 


1191 


32 


1.18 


1.80 


0.321 ±0.042 


4.64 


11.17 ± 1.46 


14.2+^ : « 



Note: Clustering results for samples with N < 500 are not shown due to poor quality statistics for our field area. 
The ERO types are labeled as E for passive EROs and SF for star-forming EROs. 
% total is the fraction of combined EROs in a particular ERO-type sub-population. 



Table 3 

Summary of combined ERO clustering results (ONLY galaxies with z phot information). 





ERO cut 


Mag. limit 


ERO type 


N 


z 


7 


w(l') 


X 2 /d-o.f 


rn (h 1 Mpc) 


M ha , (XlO 12 Mr) 

J ft tit V V \Zs s 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


ALL 


2539 


1.11 


1.89 ±0.07 


0.357 ±0.033 


3.25 


9.88 ±1.68 


9.7+ti 

-3.8 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


ALL 


7464 


1.14 


1.88 ±0.05 


0.254 ±0.019 


1.88 


8.74 ±1.21 


6.6 + " 

-2.3 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


ALL 


15380 


1.16 


1.76 ±0.06 


0.198 ±0.016 


1.56 


8.64 ±1.22 


6.3 + ?'« 

-2.4 


I- 


[3.6] > 5.0 


[3.6] < 18.0 


ALL 


22754 


1.17 


1.73 ±0.05 


0.158 ±0.012 


2.24 


7.89 ± 1.18 


4.8 + 2« 


I- 


[3.6] > 5.0 


[3.6] < 18.5 


ALL 


25703 


1.17 


1.76 ±0.05 


0.145 ±0.011 


2.80 


7.39 ±1.04 


3.9 + f ; o 


R 


-K s > 5.0 


K s < 17.9 


ALL 


2155 


1.05 


1.68 ±0.07 


0.402 ±0.043 


3.68 


14.36 ±2.65 


30.7+21 i 

— 1 j.S 


R 


-K s > 5.0 


K s < 18.4 


ALL 


5395 


1.09 


1.76 ±0.06 


0.312±0.027 


1.94 


11.67± 1.85 


16.2+" 


R 


-IQ > 5.0 


Ks < 18.9 


ALL 


10795 


1.12 


1.69 ±0.05 


0.239 ±0.021 


1.87 


10.90±1.85 


13.1+82 


R 


-K, > 5.0 


K s < 19.4 


ALL 


15209 


1.13 


1.70 ±0.06 


0.205 ±0.017 


1.37 


9.80 ±1.45 


9.4+^ 


I 


-K s > 4.0 


K s < 17.9 


ALL 


940 


1.08 


1.90±0.13 


0.290 ±0.046 


1.96 


10.29 ±2.48 


11.0+^3 


I 


-Ks > 4.0 


Ks < 18.4 


ALL 


2653 


1.15 


1.80 ±0.07 


0.353 ±0.034 


2.47 


12.12±2.22 


18.2+124 


I 


-Ks > 4.0 


Ks < 18.9 


ALL 


6133 


1.18 


1.73 ±0.06 


0.260 ± 0.023 


3.09 


10.71 ±1.83 


12.4+™ 


I 


-K s > 4.0 


Ks < 19.4 


ALL 


8988 


1.19 


1.75 ±0.06 


0.217±0.018 


1.49 


9.48 ±1.56 


8.5+5-2 


I- 


[3.6] > 5.0 


[3.6] < 16.5 


ALL 


2539 


1.11 


1.80 


0.381 ±0.030 


3.27 


11.19±0.88 


14 - 2+ 32 


I- 


[3.6] > 5.0 


[3.6] < 17.0 


ALL 


7464 


1.14 


1.80 


0.273 ±0.016 


1.91 


9.79 ±0.57 


9.4+J-8 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


ALL 


15380 


1.16 


1.80 


0.190±0.010 


1.75 


8. 16 ±0.43 


+0 8 
—0.8 


I- 


[3.6] > 5.0 


[3.6] < 18.0 


ALL 


22754 


1.17 


1.80 


0.147 ±0.007 


2.79 


7.17±0.34 


3 5+ - 6 


I- 


[3.6] > 5.0 


[3.6] < 18.5 


ALL 


25703 


1.17 


1.80 


0.139 ±0.006 


3.07 


7.00 ±0.30 


3 2+°-5 


R 


-K s > 5.0 


K s < 17.9 


ALL 


2155 


1.05 


1.80 


0.364 ±0.031 


4.32 


11.84± 1.01 


16.9+« 


R 


-Ks > 5.0 


Ks < 18.4 


ALL 


5395 


1.09 


1.80 


0.301 ±0.019 


2.12 


10.97 ±0.69 


13.4+2* 


R 


-Ks > 5.0 


K s < 18.9 


ALL 


10795 


1.12 


1.80 


0.214±0.012 


2.86 


9.21 ±0.52 


7 7+ 1 - 3 
'■'-1.3 


R 


-Ks > 5.0 


Ks < 19.4 


ALL 


15209 


1.13 


1.80 


0.185 ±0.010 


2.18 


8.45 ±0.46 


6.6+3-3 


I 


-Ks > 4.0 


Ks < 17.9 


ALL 


940 


1.08 


1.80 


0.297 ±0.048 


2.00 


11.53 ± 1.86 




I 


-Ks > 4.0 


Ks < 18.4 


ALL 


2653 


1.15 


1.80 


0.353 ±0.028 


2.47 


12.12±0.96 


18.2+4;* 


I 


-Ks > 4.0 


Ks < 18.9 


ALL 


6133 


1.18 


1.80 


0.243 ±0.016 


3.43 


9.64 ±0.63 


9.0+20 


I 


-Ks > 4.0 


Ks < 19.4 


ALL 


8988 


1.19 


1.80 


0.207 ±0.012 


1.73 


8.83 ±0.51 


6.8+j2 



Note: Clustering results for samples with N < 500 are not shown due to poor quality statistics for our field area. 
These samples contain only objects that have z p hoi information. 



Table 4 

Summary of type-separated ERO clustering results (ONLY galaxies with z p hot information). 





ERO 


Mag. 


ERO 


N 


% 


z 


7 




v 2 n f 




m halo 








type 




total 










yll ivipc J 


( Y 10 12 Mr^ 


1 — 


r 61 ~> 5 n 


R 61 ^ 1 6 S 


E 


1849 


73 


1.11 


i .70 zn u.uo 


500 4- 045 


2.95 


i 34 + 1 79 
i u.jt zlz i . i y 


1 1 1+ 71 


J— 


r fii > 5 o 

|_*7.UJ ^ .YJ 


R 61 < 1 6 5 


SF 


690 


27 


1.09 




758 + 063 


6.83 


11 os -|_ 3 20 


17 1+I8.6 
U - L -1Q.5 


I 


R 61 >> S 


r 6i <" 1 7 n 

[J .UJ \ 1 / .YJ 


E 


5003 


67 


1 14 


9 01 -4- O 06 
Z..U1 zrz u.uu 




2.55 


Q « -L- 1 11 
o. JVJ 1.11 


6 1+ 27 


I 


r 6i ~~> s n 


r 6i <? 17 n 


SF 


2461 


33 


1 14 


i _|_ n no 

1 .UJ ZTZ U.U7 


U. 1 OO ZC U.UZ.U 


5 68 


9 95 + 2 17 


9 9+83 


I- 


[3.6] > 5.0 


[3.6] < 17.5 


E 


9697 


63 


1.16 


1.90 ±0.05 


0.249 ±0.01 8 


1.28 


8.30+ 1.14 


5 6 +2 | 


I- 


R 61 > 5 


R 61 < 1 7 5 


SF 


5683 


37 


1.16 


1.88 ±0.08 


0.142 ±0.014 


7.01 


6.95 ± 1.33 


-i 2+2.4 
■'■''-1.5 




-K s > 5.0 


K s < 17 9 


E 


1761 


82 


1.04 


i q? + f) 08 

1 . 7 £. _1_ U . \J O 


401 + o 039 


3.49 


1 104 ± 1 .94 


13 6 +8 ' ) 

—5 9 


R 


-K, > 5.0 


if, < 17.9 


SF 


394 


18 


1.06 












R 


-K s > 5.0 


7f, < 18.4 


E 


3804 


71 


1.08 


i so + () ()7 

± .YJ y _1_ U . V./ / 


327 + 028 


3.15 


10 35 + | 56 


1 1 2 +5<) 


R 


-K s > 5.0 


K s < 18.4 


SF 


1591 


29 


1.11 


1.61 ±0.08 


0.306 ± 0.042 


2.36 


13.71 ±3.09 


26 7+ 23 - 2 
'"■'-13.6 


J 


-K, > 4.0 


if, < 17.9 


E 


802 


85 


1.07 


2.11 ±0.15 


318 + 051 


1.51 
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Note: Clustering results for samples with N < 500 are not shown due to poor quality statistics for our field area. 
The ERO types are labeled as E for passive EROs and SF for star-forming EROs. 
% total is the fraction of combined EROs in a particular ERO-type sub-population. 
These samples contain only objects that have z p ha information. 



